Effects of finite temperature on the Mott insulator state 



O 

o 

(N 

c 

a 

m 

(N 



5-1 

43 



c 

o 
o 



> 
in 
i> 
o 

o 

o 

-I— > 



S3 

o 
o 



X 



Guido Pupillo 1 ' 2 ' 3 , Carl J. Williams 1 , and Nikolay V. Prokof'ev 4 5 

1 National Institute of Standards and Technology, Gaithersburg, MD 20899 
2 Department of Physics, University of Maryland, College Park, MD 20742 
3 Institute for Quantum Optics and Quantum Information, 6020 Innsbruck, Austria 
4 Department of Physics, University of Massachusetts, Amherst, MA 01003 
5 Russian Research Center "Kurchatov Institute", 123182, Moscow 

We investigate the effects of finite temperature on ultracold Bose atoms confined in an optical lattice plus 
a parabolic potential in the Mott insulator state. In particular, we analyze the temperature dependence of the 
density distribution of atomic pairs in the lattice, by means of exact Monte-Carlo simulations. We introduce 
a simple model that quantitatively accounts for the computed pair density distributions at low enough temper- 
atures. We suggest that the temperature dependence of the atomic pair statistics may be used to estimate the 
system's temperature at energies of the order of the atoms' interaction energy. 



I. INTRODUCTION 



Great progress has been achieved in the coherent control 
and manipulation of ultracold atoms and molecules. In a se- 
ries of recent experiments, several groups were able to demon- 
strate the loading of an atomic Bose-Einstein condensate into 
the lowest vibrational level of single wells of an optical lat- 

En n n r~i 
,I2L|3|,|J.|3|]- In a remarkable experiment, M. Greiner et 

ai, demonstrated the capability of inducing a reversible 
quantum phase transition between a superfluid and an insu- 
lator state for bosonic atoms by varying the intensity of the 
trapping laser beams and therefore the depth of single poten- 
tial wells. In a series of recent experiments, the same group 
demonstrated the capability of creating an array of quasi one 
dimensional tubes and induce an effective one dimensional 
superfluid/insulator transition, with on-site peak density equal 
to one olIOI- Several proposals suggest this system has 
potentially fundamental applications in the field of quantum 
computation, as the realized array of atoms in the insulat- 
ing re gime can be used to create an ideal register for qubits 

A key observable in experiments with ultracold atoms is 
the interference pattern observed after a certain time of flight 
following the release of atoms from all trapping potentials. 
The interference pattern is proportional to the distribution of 
atomic momenta 114 ll 5l llq [1711 . and provides direct infor- 
mation on the amount of phase coherence present in the sys- 
tem. Information extracted from the interference pattern is 
routinely used to measure the system temperature in the su- 
perfluid regime, but cannot be easily used in the strongly inter- 
acting regime, where both finite temperature and interactions 
compete to fill up the single-particle Brillouin zone. In fact, 
in time-of-flight images the measurement precision is on the 
order of the level spacing to the second lattice band, while the 
relevant energy scales in the strongly interacting regime are 
the onsite interaction energy and the hopping amplitude. The 
latter are typically at least one order of magnitude smaller than 
the lattice level spacing fisll . Thus, apart of a scheme based 
on fermionization of the many -body wavefunction |6], to date 
there are no available schemes for the estimation of temper- 
ature at energies of the order of the interaction energy in the 
strongly correlated regime. 



Recently, time-of-flight images have provided evidence for 
the mixing of particle-hole excitations in the ground state in 
the Mott regime Ref. 0,0]. This mixing has important con- 
sequences for the use of the Mott state as a quantum register. 
In fact, the presence of empty and doubly occupied sites in- 
troduces unwanted errors in the register initialization, which 
have to be corrected for at the cost of overheads in compu- 
tational resources and gate times. A characterization of the 
finite-temperature population of particle-hole states in realis- 
tic experimental conditions is thus desirable. 

In this work we analyze the effects of finite temperature 
on the Mott insulator state in the presence of an external 
parabolic potential. We focus on the effects of temperature 
on the density distribution of doubly occupied sites, showing 
that at finite temperature the latter has a characteristic site- 
dependent Gaussian profile due to the presence of the external 
potential. We derive the functional form of this temperature 
dependence by introducing an effective single-particle model 
for the many-body spectrum at energies of the order of the in- 
teraction energy, and we compare our model's predictions for 
the atomic pair density distribution to exact quantum Monte- 
Carlo results for realistic experimental parameters of the trap- 
ping potentials, finding good agreement. We suggest that the 
temperature dependence of the density distribution of doubly 
occupied sites may be used to estimate the system's tempera- 
ture at energies of the order of the interaction energy. 

The presentation of the results is organized as follows. In 
Sec.|n|we review the physics of the Mott insulator state and 
the low-energy spectrum of the Bose-Hubbard Hamiltonian 
both in homogeneous and inhomogeneous lattices. In Sec.lIIII 
we introduce an effective model for the many-body spectrum 
at energies of the order of the interaction energy, and provide 
an analytical expression for the temperature dependence of 
the density distribution of atomic pairs at low enough temper- 
atures. In Sec. II VI we validate our model by comparing the 
analytical results to exact quantum Monte-Carlo simulations. 
In Sec. [V] \ve address the stability of the Mott insulator state 
against finite temperature population of particle-hole states. In 
Sec. lVll we discuss the possibility of estimating the system's 
temperature by collecting statistics of the presence of atomic 
pairs in the lattice. Finally, the conclusions are presented in 
Sec. ITU 
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II. THE MOTT INSULATOR STATE 

The Bose-Hubbard Hamiltonian describes the system's dy- 
namics when the lattice is loaded such that only the lowest 
vibrational level of each lattice site is occupied 1191 



N = M. In the presence of the trap these states are not de- 
generate, spanning an energy range which may be larger than 
U. We redefine as \T) = n^^^jlO) the lowest en- 
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Here ctj is the bosonic annihilation operator of a particle 
at site j and nj = OjCij. The site-dependent e(j) mod- 
els an external magnetic quadratic potential, while J and 
U are the tunneling and on-site interaction energies respec- 
tively. For deep lattices with potential energy V(x) = 
Vcos 2 (kx), the tunneling energy is approximated by J = 

A/{^)E R {V/E R f/ i e- 2 ^ v / ER , where E R = (hk) 2 /2m 
is the recoil energy, k the light wave vector, and m the atomic 
mass. The on-site interaction energy is due to ground state 
collisions between atoms each in the motional state 4>(x) and 
is given by U 
tering length and 



J dx.\(j>(x)\ 4 , with a s the s-wave scat- 
x) a Wannier state. J and U are assumed 
to be site independent | ljj . 

For a homogeneous one dimensional lattice (e(j) — 0, Vj) 
of M sites with periodic boundary conditions at zero temper- 
ature T, an insulating state, the Mott state, occurs only if the 
number of particles N equals M and the interaction energy 
dominates over the tunneling energy. If y/~N J/U -C 1 the 
ground state has approximately energy E g — —ANJ 2 /U and 
is given by 



\Hf g ) = a(\T)+2VNJ/U\S)), 



(2) 



where the normalization constant is c 
AN {J/U) 2 )- 1 / 2 . The Fock state \T) = U 
one particle per site and zero energy, and the symmetrized 



ergy state and set its energy to zero 1 20] . The state |T) is 
the ground state of the system to zeroth order in J, while the 
true ground state has coherences due to tunneling of a hole 
in either one of sites ±(JV — 1) /2 and mixing of Fock states 
\Sn l7n ) = a^amlT) / y/2 with an atomic pair and a hole at sites 
n and m respectively, with n,m < \(N — 1)/2|, in analogy 
to the homogeneous case. States \S ntm ) are not degenerate, 
due to the trap presence. States \S 7um ) with m = n ± 1 dom- 
inate the mixing into the ground state, and for £1 <C U the 
amplitude of mixing approaches the homogeneous system's 
one V2aJ/U. 

In general, the |T)-state is coupled via the BH- 
Hamiltonian to |S nim )-states with a matrix element of order 
( J/?/)'™ - " 1 ', which is negligible for \n — m\ large enough in 
the insulating regime. On the other hand, a temperature of the 
order of the energy difference between |T) and \S n;m ) -states 
allows for population of the latter. In particular, the lowest 
energy Fock state \S n ,m) nas two particles in the central site 
of the lattice and a hole at site \(N — l)/2\. Its energy is 
A = U — e((N — l)/2). This state is an eigenstate of the 
£?i?-Hamiltonian for J = only. For any finite J, eigen- 
states are a superposition of many different Fock states. This 
superposition lowers the energy A of an amount W which is 
of the order of a few J. As a consequence, the energy A of the 
lowest-energy eigenstate with two atoms in a site is reduced 
to approximately A = A — W. In the following we introduce 
a simple model for the spectrum of the many-body system at 
energies of order A, which we find to be the relevant energy 
scale for the double occupation of a site in the lattice fllll . 
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state \S) = l/V2NY,f=i( a ]a J+1 + a) +x aj)\T) has energy III. EFFECTIVE SINGLE-PARTICLE SPECTRUM 



U. Indeed, in the homogeneous system all Fock states with an 
empty site and an atomic pair in another site are degenerate 
with energy U. The latter is roughly the lowest excita- 
tion energy for the many-body system in the homogenous 
commensurately filled case | lJI. 

The probability of creating exactly N = M particles in an 
experiment is negligible. An insulator state can be recovered 
if an external quadratic potential is present. Then, e(j) = Qj 2 
where il = m(ir/k) 2 Lj 2 -,/2 is a characteristic energy scale for 
a trap of frequency lot- At T = 0, a Mott state with average 
occupation one in the center of the trap can be realized for 
VNJ/U < 1 if 



N < M, U > e((N - l)/2), 



nN > J. 



(3) 



The last two inequalities insure that multiple particle occu- 
pancy and tunneling of holes from the borders to the cen- 
ter of the lattice are suppressed, respectively 1121 . The first 
inequality in Eq. [3] simply states the number of wells be 
larger than the number of atoms. As a consequence, there 
are Ml/((M — N)\Nl) Fock states with maximal on-site oc- 
cupation one, at variance with the single state of the case 



Neglecting for a moment mixing of |5 n>m )-states in low ly- 
ing modes, we expect that the density distribution of atomic 
pairs as a function of the position in the lattice at finite tem- 
perature T has a a Gaussian profile for the following argu- 
ment. Given a lattice with M wells, N = M + 1 particles 
and U 7$> J, the ground state has roughly one particle per 
site and an extra particle at the center of the trap, forming 
an atomic pair. For weak enough quadratic traps, £1 <C J, 
tunneling of the extra particle is not suppressed over a cer- 
tain number of lattice sites at the center of the trap where the 
external potential is essentially flat, meaning that the atomic 
pair acts approximately as a conventional harmonic oscilla- 
tor with effective mass = 1/(4 J) and trapping frequency 
lu* = 2a/ 217 J. Here the lattice constant a — Txjk and 
h = h/2ir have been set to unity, where h is the Planck con- 
stant. The tight binding spectrum is therefore approximated 
by an harmonic oscillator spectrum whose level spacing is 
with an error of order 0(il/J). If T > ui^/kg, it can 
be shown that the quantum density distribution P(j) for the 
atomic pair becomes equivalent to the classical distribution 
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which is proportional to 



0.01 



-pnf 



(4) 



where (3 — l/(fcsT), and fc# is Boltzmann constant. The 
distribution P(j) as a function of lattice site j, Eq. |4] is 
therefore a Gaussian whose width is xq = y/l/(2(3tt). 

We postulate that for realistic situations where M > N 
harmonic oscillator states actually approximate the incom- 
plete basis set defined on the central N sites of the lattice in 
the energy range A < E < 3U - 2e((N - l)/2), where 
3U — 2e((A — l)/2) is a characteristic energy for population 
of states with three particles in one of the central N — 2 sites. 
Apart of a normalization factor, we therefore expect that the 
density distribution as a function of the position j in the lattice 
has the same Gaussian profile as in Eq.|4]with an exponential 
suppression due to the energy shift A 



P(j) oc e 



(5) 



The width xo of the density distribution can be there- 
fore directly related to the temperature of the sample, for 
ksT < A. In fact, while for a temperature of order A states 
are populated with more than one pair in the central sites, 
the various pairs approximately behave like non-overlapping 
harmonic oscillators, meaning that the probability of finding 
two extra particles in the same site is strongly suppressed. 
Therefore, the density distribution at a single site should 
still approximately maintain a Gaussian profile of width xq. 
Moreover, the height of the Gaussian peak is an indication of 
A, and therefore of the number of particles N, For vanishing 
values of J, A is approximately A, and the peak's height is 
easily computed. For J > 0, the shift W = A — A is of 
the order of half of the pair's band-width. Because the latter 
is 8J, W is approximately 4J. Because the peak's height 
also depends on the normalization, we expect Eq. [5] to well 
approximate the value of the density peak for fc^T <C A, 
where the normalization is almost one. 




FIG. 1: Density distribution of atomic pairs P as a function of the 
lattice index j. Continuous and dashed lines are numerical results 
for T = 3J/fcs and T = 5J/fc_g respectively. Dotted lines are 
analytical curves for the same temperatures. The arrow indicates the 
zero-temperature mixing of | S n ,n±i) -states. HereiV = 101, [//J = 
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The modeling of the excited states of the many-body sys- 
tem with an effective harmonic oscillator spectrum for atomic 
pairs is one of the central results of this paper. In the follow- 
ing we verify the accuracy of this picture by comparing the 
model's predictions for the atomic pair density distribution to 
exact Monte-Carlo results. 



IV. NUMERICAL SIMULATIONS 

We perform numerical simulations of the density distribu- 
tion of atomic pairs in the lattice as a function of temperature 
for some realistic experimental parameters. We employ a 
quantum Monte-Carlo code based on the continuous-time 
Worm algorithm |21]. A sample of A = 101 atoms of 87 Rb is 
trapped in a lattice with wavelength A = 785 nm and parallel 
and transverse confinements Vj| = 15Er and V± = 40.5i?R 
respectively. Using a s — 5.6 nm, the interaction energy in 
frequency units is U/h = 3.340 kHz and U/J = 120. The 



FIG. 2: Density distribution of atomic pairs P as a function of lattice 
index j, for T — 12J/fcs. Dashed and dotted lines are numerical 
and analytical results respectively. Here N = 101, U / J — 120, 
Q/J = 0.0374226, and A/J = 24. 



external magnetic trapping frequency is lot ~ 27r x 40 Hz and 
therefore the number of atoms N fixes e(N — l)/2 w 0.8U, 
and thus A « 24 J. We have chosen these parameters because 
they are experimentally feasible and satisfy Eq. [5] so that 
at zero temperature a Mott state is formed with one particle 
occupying each one of the central N sites of the lattice. In par- 
ticular, zero temperature mixing of basis states with an empty 
site in one of the central N sites of the lattice is suppressed, 
the largest amplitude of mixing corresponding to the outer 
most trapped particle tunneling to the adjacent empty site. 
Calculations have been performed for M 3> N, so that atoms 
never reach the border of the lattice. We note that N = 101 
compares well with the typical number of atoms present in 
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FIG. 3: (Color online) Width xo of the density distribution of atomic 
pairs as a function of fesT/A. The width xo is in units of the lattice 
constant a. Circles(red) are numerical values, while the continuous 
line is the analytic curve (k B T /2Q,) 1/2 . Here N = 101, U/J = 
120, fi/J = 0.0374226, and A/ J = 24. 
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the one-dimensional tubes of Refs. J3 H B H H II H For 
example, N is equal to 80 in the central tube of Ref. 01, 
which is the largest number of atoms for the entire arrayof 
tubes, while N is as small as 20 in the central tube of Ref. 1 6] . 

In Fig.[2the atomic pair density distribution P is plotted as 
a function of the site index j for Tks/J — 3 and 5. Con- 
tinuous and dashed lines are numerical results, while dotted 
lines are analytical curves. For Tfce/J = 3, lower curves, 
the numerical solution shows essentially a flat density distri- 
bution throughout the central N sites, with a shallow gaus- 
sian peak at the center. The flat distribution corresponds to 
the zero temperature residual mixing of |S' n! „±i)-states into 
the ground state, characteristics of the Mott state. The results 
of Ref. [7] are consistent with the observation of this mix- 
ing. As f2 <C U, corrections to the density distribution due to 
the external trapping potential are not distinguishable on the 
scale of the graph, and the on-site density matches the homo- 
geneous system's value 2 x 2(J/U) 2 w 0.000278. The factor 
of two in front of last expression is due to the fact that the 
extra atom can tunnel from the left or right. The latter con- 
stant has been added to the analytic curves, and is indicated 
by an arrow in Fig.^ Finite atomic pair density due to zero- 
temperature mixing of |5 rii „±i)-states is a direct signature of 
the creation of a Mott state. For sufficiently low temperatures, 
selective probe of the density distribution in lattice sites j with 
j ^s> xq can give direct in situ evidence of the formation of 
the Mott plateau in the center of the trap, while the shallow 
peak around j = measures finite temperature population of 
1 5 n , m ) -states. 

For Tks I J = 5 the Gaussian peak is more evident on the 
scale of the graph. We observe that numerical and analytical 
curves nearly overlap. In particular, the widths of the Gaus- 
sians perfectly match, while the height of the peak is slightly 
underestimated by the analytic curve. This is due to the fact 
that A is only approximately 24 J for the chosen parameters, 



FIG. 4: Local densities and fluctuations as a function of the site in- 
dex j for k B T/J = 3(a), 5(6), 12(c), and 24(d). Here, A/J = 24. 
Continuous and dotted lines are the numerical densities and fluctua- 
tions, respectively. 



and the difference is amplified by the exponential function. 
A larger disagreement in the peak height is observable for 
ksT/J= 12 in Fig. [2] where the analytic curve largely over- 
estimates the exact result. This is expected, as for higher tem- 
peratures the normalization constant differs more and more 
from one, therefore suppressing the peak height in the exact 
solution. 

In Fig. |3] the width xo is plotted against the temperature, 
up to T = A/ks- Circles are numerical results, while the 
continuous line is the analytic solution yj ksT '/ (20). The 
plot shows an almost perfect agreement between numerical 
and analytical results for fc^T < A/4. Deviation for higher 
temperatures is due to population of states with a large 
number of atomic pairs. For fc^T = A the deviation is still 
less than a lattice constant. We consider the good agreement 
between the analytic model and the numerical results over the 
entire range ksT < A as a strong indication of the validity 
of the harmonic oscillator model. 



V. FINITE-T DEPLETION OF THE MOTT STATE 

A large population of states with two particles in a site will 
ultimately totally deplete the Mott state, thus making any ap- 
plication of the Mott insulator state to quantum computation 
schemes impossible. In the remainder of this section we inves- 
tigate the role of the energy scale A in the thermal depletion of 
the Mott insulator state, by numerically analyzing the atomic 
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density profile and onsite number fluctuations. 



In Fig.[4]the on-site atomic density 
and the number fluctuations (An/ 



j) (continuous lines) 



(dot- 



n?) 



ted lines) are plotted as a function of the site index j, for 
k B T/J = 3(a), 5(6), 12(c), and 24(d). For k B T/J = 3 
the density profile is flat, which corresponds to an extended 
Mott state around j = 0. The number fluctuations are also 
flat throughout the lattice, except for two large peaks around 
±N/2 and a small bump at the trap center. The flat regions 
correspond to the zero-temperature mixing of particle-hole 
states l22ll . which is approximately given by 2V2J/U, 
while the large peaks are due to the tunneling of particles to 
unoccupied sites, and are equal to flN II 6l 1221 . The small 
bump at the trap center corresponds to the finite-temperature 
population of high-energy states with two particles in a site, 
analogous to the discussion of Fig. ^ Panels (b) and (c) show 
that for higher temperatures atoms tend to accumulate at the 
trap center, therefore decreasing the extension of the region of 
constant density. For k B T = 24 J = A the Mott state has 
completely disappeared. 

In Ref. 12311 the effects of finite temperature on the stabil- 
ity of the Mott insulator state have been studied by means of 
classical Monte-Carlo simulations, finding that well defined 
Mott insulating plateaux are present for temperatures at most 
of the order of UJ (lQk B ). The discussion above shows that, 
for any fixed U and Q, this is true only in the limit of a small 
number of particles TV. In fact, we showed that the rele- 
vant energy scale for population of doubly occupied site is 
A = U - Q((N - l)/2) 2 , which is iV-dependent. In par- 
ticular, for a system with TV = 110 instead of TV = 101, 
we have A w 9 J. In this case, a temperature of the order 
of XJ j (lOfee) = 12 J/k B would be larger than A and would 
therefore completely deplete the Mott insulator state. 



VI. TEMPERATURE ESTIMATION 

As mentioned above, it is notoriously difficult to measure 
temperatures of the order of the interaction energy in the 
strongly interacting limit. 

On the other hand, we have shown above that the atomic 
pair density distribution has a clear temperature dependence. 
Here we suggest that in principle this temperature dependence 
may be used to estimate the system's temperature at energies 
of the order of A, by devising a probe that is sensitive to the 
presence of multiply occupied sites in the lattice. For this 
purpose, two different kinds of experiments may be consid- 
ered. The first possibility is to design an experiment which 
is solely sensitive to the total number of pairs in the lattice. 
The latter can be calculated numerically as a function of tem- 
perature by means of quantum Monte-Carlo simulations, and 
the results of the simulations can then be compared to the ex- 
perimental data to extract a temperature estimation. An ex- 
ample of these calculations for the trapping parameters given 
above and N = N to t — 101 is given in Fig. |5] where the 
total number of pairs N2 is plotted as a function of the tem- 
perature. In the plot, the circles (red) are the quantum Monte- 



Carlo results for Ntot = 101, while the squares (blue) and 
the diamonds (green) correspond to a total number of parti- 
cles Ntot — 95, and 107, respectively. Each curve is nor- 
malized to N = 101. The curves for Ntot — 95, 107 are 
shown in order to take into account a plausible experimental 
uncertainty on the total number of atoms Ntot' In fact, it is 
experimentally challenging to control the total number of par- 
ticles Ntot to better than 5 percent, and the uncertainty on the 
total number of particles is likely to dominate the experimen- 
tal errors Il24ll . Then, assuming unit efficiency of atomic pair 
detection, for each measured N2 the temperature is easily es- 
timated. For example, for the system considered above, where 
Ntot — 101, if 5 pairs (N%/N w 0.04) are measured exper- 
imentally, a simple analysis of Fig.|5]yields a most-probable 
temperature of about 14.4 J/k B (0.6 A/k B ). Thus, by tak- 
ing into account the uncertainty on the total number of atoms 
(black-dashed line in the plot), the final estimate for the tem- 
perature is T = (14.4 ± 4.0) J/k B . 

We note that the discussion above refers to the determina- 
tion of temperature in a single experiment. If one wants to 
determine average system's temperatures, the measurement 
can be repeated many times, and the error bars are obtained 
by assuming a shot-to-shot Poissonian fluctuation of N to t, in 
analogy to the discussion above. However, in the remainder 
of this section we discuss an experiment that allows for the 
estimation of average system's temperatures and it is largely 
independent of the fluctuations of N to t (although it is admit- 
tedly difficult to implement). 

In principle various schemes may be used to experimen- 
tally detect the presence of pairs in the lattice. Here, we 
propose to collect statistics of atomic pairs through resonant 
photo-association of two atoms in a single well to a molecular 
excited state, followed by ionization of the formed molecules 
and detection of the emitted ions. In fact the latter can 
be performed with almost unit efficiency ll24ll . For this 
experiment, the molecular excited state should be chosen 
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FIG. 5: (Color online) Number of atomic pairs N2 as a function of 
kBT/A. Here, N2 is in units of N = 101. The squares(blue), cir- 
cles(red), and diamonds(green) correspond to systems with Ntot = 
95, 101, and 107 atoms, respectively. For all curves, A/J = 24. 
The dashed line is the error bar for the temperature estimation, given 
a measured number of pairs equal to 5 (N2/N ~ 0.04). The dotted 
lines are a guide to the eye. 
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to be far from atomic dissociation, meaning that the atomic 
probability of spontaneous emission is not enhanced by the 
photo-associative laser. 

The second kind of experiment that may be envisioned is 
one where the probe has a high spatial resolution, so that it 
is possible to determine the position of the atomic pair in the 
lattice. In this case, information on average system's tempera- 
tures may be extracted by measuring the width ieo, which has 
been shown above to depend only on the temperature and the 
known trap geometry. The advantage of this experiment is that 
xo is insensitive to N to t- Thus, calibration of average system's 
temperatures may be performed by accumulating statistics of 
pair detection on successive experiments with the same trap- 
ping potentials, without worrying about possible fluctuations 
of Ntot from one experiment to the next. 

Although a careful discussion of the sources of uncertainty 
on the measurements is out of the scope of the present 
analysis, we expect that a significant uncertainty on the tem- 
perature measurement is likely to come from the finite spatial 
resolution of the probe-field. In fact, for typical experimental 
setups it is not possible to focus a photo-associative laser 
onto a single site. The intensity of the photo-associative laser 
has a Gaussian profile whose width <j£ is typically of the 
order of the light's wavelength or larger. This means that 
the laser intensity may not be negligible on a few lattice 
sites, depending on the ratio between the wavelength of 
the photo-associative laser and the lattice constant a. A 
convolution between the pair density-distribution Gaussian 
and the laser-intensity Gaussian provides an estimate of the 
experimentally measured width x' Q , which is thus given by 

x 'o = V x o + °l ■ We note that the difference between xo 
and x' decreases with increasing temperature. For example, 
for a reasonable value ul = 3a, x' — xo is approximately 
0.6 a for ksT — 0.2A, while it drops to x' — xq rs 0.3a for 
fcgT = 0.6A. 



Finally, we note that it is straightforward to generalize the 
calculations above for the array of one-dimensional tubes of 
Refs. lyl |3l |4|, |5l |a [3, |8|] . Assuming an initial Thomas-Fermi 
distribution of the atoms in the three dimensional system, the 
calculations above may be repeated for each tube separately 
and the results combined (see also Ref. |6]). 



VII. CONCLUSIONS 

In summary, we have studied the effects of temperature on 
the Mott insulator state, with an emphasis on the finite tem- 
perature population of states with two atoms in a site. We 
introduced a model for the energy spectrum of the trapped 
many-body system in the strongly correlated regime, which 
allowed us to derive the functional form of the dependence of 
the density of atomic pairs upon temperature. We discussed 
the depletion of the Mott state due to the population of states 
with multiple on-site occupancy. Finally, we suggested that 
the temperature dependence of the density distribution of dou- 
bly occupied sites may be used to estimate the system's tem- 
perature at energies of the order of the interaction energy. 

Note added: After this work was completed, the results 
of an experiment analogous to one of those proposed here 
have been presented in Ref. ll25ll . where the temperature of 
a fermionic Mott state has been estimated. 
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